arXiv:1501.07658v2 [cs.IT] 18 May 2015 


1 


Robust Transceiver Design for MISO Interference 
Channel with Energy Harvesting 

Ming-Min Zhao, Yunlong Cai, Qingjiang Shi, Benoit Champagne, and Min-Jian Zhao 


Abstract —In this paper, we consider multiuser multiple-input 
single-output (MISO) interference channel where the received 
signal is divided into two parts for information decoding and 
energy harvesting (EH), respectively. The transmit heamforming 
vectors and receive power splitting (PS) ratios are jointly designed 
in order to minimize the total transmission power subject 
to both signal-to-interference-plus-nolse ratio (SINK) and EH 
constraints. Most joint heamforming and power splitting (JBPS) 
designs assume that perfect channel state information (CSI) is 
available; however CSI errors are inevitable in practice. To over¬ 
come this limitation, we study the robust JBPS design problem 
assuming a norm-bounded error (NBE) model for the CSI. Three 
different solution approaches are proposed for the robust JBPS 
problem, each one leading to a different computational algorithm. 
Firstly, an efficient semidefinite relaxation (SDR)-based approach 
is presented to solve the highly non-convex JBPS problem, where 
the latter can be formulated as a semidefinite programming (SDP) 
problem. A rank-one recovery method is provided to recover 
a robust feasible solution to the original problem. Secondly, 
based on second order cone programming (SOCP) relaxation, we 
propose a low complexity approach with the aid of a closed-form 
robust solution recovery method. Thirdly, a new iterative method 
is also provided which can achieve near-optimal performance 
when the SDR-based algorithm results in a higher-rank solution. 
We prove that this iterative algorithm monotonically converges 
to a Karush-Kuhn-Tucker (KKT) solution of the robust JBPS 
problem. Finally, simulation results are presented to validate the 
robustness and efficiency of the proposed algorithms. 

Index Terms —MISO interference channel, beamforming, 
power splitting, semidefinite programming, second-order cone 
programming, concave-convex procedure. 

I. Introduction 

Recently, energy harvesting (EH) from the environment has 
attracted considerable interest since it offers a promising solu¬ 
tion to provide cost-effective and perpetual power supplies for 
wireless networks Q-ID- The unified study of simultaneous 
wireless information and power transfer (SWIFT) has therefore 
drawn significant attention lately, as it opens new challenges 
and possibilities in the analysis and design of transmission 
schemes and protocols. 
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The fundamental concept of SWIFT was first proposed in 
Q, which characterizes the rate-energy (R-E) tradeoff in a 
discrete memoryless channel. The study of R-E tradeoff was 
later extended to frequency selective channels multiple 
access and multi-hop channels |j^, and two-way channels 
Q- However, the above works all assume that the receiver 
can decode information and harvest energy from the same 
signal, which is impossible with existing technology. In Q, 
the authors proposed a practical receiver structure for the first 
time, and considered the R-E region and optimal transmission 
scheme of a MIMO broadcasting channel. Two practical signal 
separation schemes, namely time switching (TS) and power 
splitting (FS) were also considered. For the TS scheme, the 
transmitter divides the transmission block into two orthogonal 
time slots, one for transferring power and the other for 
transmitting data. For the FS scheme, the received signal is 
split with an adjustable FS ratio, where the stream with power 
ratio 0 < p < 1 is used for information decoding (ID) and the 
other stream with power ratio 1 — p is used for EH. 

The works in and 0 considered beamforming de¬ 
signs with separate information/energy receivers. Specially, 
@ studied the robust beamforming problem for the multi¬ 
antenna wireless broadcasting system with SWIFT, under the 
assumption of imperfect channel state information (CSI) at 
the transmitter. In |]^ the authors investigated the optimal 
beamforming strategy to maximize the weighted sum-power 
transferred to all EH receivers subject to given minimum 
signal-to-interference-plus-noise ratio (SINR) constraints at 
different ID receivers. The work in ||^ derived the optimal 
power splitting rule at the receiver to achieve various trade¬ 
offs between the maximum ergodic capacity for information 
transfer and the maximum average EH for power transfer. In 
various practical receiver architectures for SWIFT were 
investigated, where a new integrated information and energy 
receiver design was proposed. The use of SWIFT for relay 
systems was considered in pO|-p^. 

The work d) considers a multiuser MISO downlink 
system with SWIFT, where the total transmission power at 
the base station (BS) is minimized subject to given SINR 
and EH constraints. An optimal solution is proposed based 
on semidefinite relaxation (SDR) along with low-complexity 
suboptimal designs. Some recent studies have focused on 
SWIFT in the context of multi-antenna interference channels 
ig-iiizl- The work IB investigates optimal transmission 
strategies and mode scheduling methods for a two-user MIMO 
interference channel with EH, while extends this study 
to the AT-user MIMO interference channel. A number of 
recent studies focus on the FS approach and consider the joint 
beamforming and power splitting (JBFS) design problem in 
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MISO interference channel, where the downlink receivers are 
characterized by both SINR and EH constraints. Compared 
with the conventional beamforming design, JBPS design is 
much more challenging due to the coupling between the 
beamforming vectors and PS ratios. Hence, the corresponding 
research area has become quite active and several algorithms 
have been recently proposed to address this problem. In (16| 
and flTt , the JBPS design is studied for a itT-user MISO inter¬ 
ference channel with the same design criterion as that in |T3). 
Specially, the work uses SDR to tackle the non-convex 
JBPS problem and prove that the SDR is tight when K — 2 
or 3. Also, various suboptimal but low complexity solutions 
based on fixed beamforming schemes and a hybrid scheme 
are provided. In flT) , the JBPS problem is reformulated as a 
SOCP problem based on an alternative method named SOCP 
relaxation and two sufficient conditions are given under which 
the relaxation is tight. Also, a primal-decomposition based 
distributed algorithm is proposed for the JBPS problem. 


In these existing works on JBPS, the CSI is assumed to 
be perfectly known. In practice however, the CSI is prone to 
errors owing to various factors, which may limit the system 
performance drastically. Hence, it is of interest to develop 
JBPS designs that are robust to CSI errors. In this paper, 
assuming norm-bounded CSI error (NBE) models for the CSI, 
we propose three new computational algorithms for solving the 
robust JBPS design problem in a AT-user MISO interference 
channel with multi-antenna transmit beamformers and single¬ 
antenna PS receivers. In the first approach, we show that the 
robust JBPS problem can be relaxed as a SDP problem based 
on SDR. A rank-one recovery method is provided to recover 
a robust feasible solution to the original problem if a high- 
rank solution is returned. In the second approach, we propose 
to formulate the original problem as a SOCP problem based 
on SOCP relaxation and the Cauchy-Schwarz inequality for 
the purpose of complexity reduction. Since the solution to the 
SOCP problem is not necessarily robust, a closed-form robust 
solution recovery method is provided. Einally, as our third 
approach, we propose a new iteration algorithm based on the 
concave-convex procedure (CCCP) which can provide near- 
optimal performance when higher-rank solutions are returned 
by the SDR-based algorithm. The convergence of this iterative 
algorithm is studied in detail and we prove that any limit point 
of the iterative algorithm is a Karush-Kuhn-Tucker (KKT) 
solution to the robust JBPS problem. Einally, simulation results 
are presented to validate the robustness and efficiency of the 
proposed algorithms. In particular, we show that the proposed 
robust algorithms can provide near-optimal performance. 


The reminder of this paper is organized as follows. We 
present the system model of the AT-user MISO interference 
channel, the channel error model and the robust JBPS problem 
in Section]^ In Section]^ the robust JBPS design approaches 
based on SDR and SOCP relaxation are developed along 
with their rank-one recovery methods and final computational 
algorithms. In Section IV we develop the CCCP-based iter¬ 
ative design algorithm for the JBPS problem and discuss its 
convergence and initialization. A detailed complexity analysis 
of the proposed algorithms is provided in Section |V] Einally, in 
Section VI computer simulations are used to compare the ro¬ 


bust performance of the proposed JBPS designs. Conclusions 
are drawn in Section Ivnl 

Notations: Scalars, vectors and matrices are respectively 
denoted by lower case, boldface lower case and boldface upper 
case letters. For a square matrix A, Tr(A), rank(A), A^, 
conj(A), A^, A~^ and A^^ denote its trace, rank, transpose, 
conjugate, conjugate transpose, inverse and pseudo-inverse, 
respectively, while A ^ 0 means that A is a positive semidef- 
inite matrix. The operator diagjSi,..., Sm} denotes a block- 
diagonal square matrix with Si,..., Sm denoting the diagonal 
square matrices. The operator vec(-) stacks the elements of a 
matrix in one long column vector, invp(a;) denotes the inverse 
of the positive portion, i.e. . || • |j denotes the Euclidean 

norm of a complex vector and | • | denotes the absolute value 
of a complex scalar. Finally, denotes the space 

of TO X n complex (real) matrices, and IR+ (®-) denotes the 
space of positive (negative) real numbers. 

H. System Model And Problem Formulation 
A. Proposed System Model 

We consider the A-user MISO interference channel where 
each transmitter, indexed by fc € /C = {1,..., A}, is equipped 
with Afc antennas and each receiver is equipped with a single 
antenna. The A transmitters are assumed to operate over a 
common frequency band and each communicates with its 
corresponding receiver using transmit beamforming. Different 
from conventional interference channels, we here consider 
PS-based receivers. The received signal at each receiver is 
split into two separate signal streams with different power 
levels, one sent to the EH receiver and the other to the ID 
receiver ||^. The system model is illustrated in Fig. [T] We 
assume that transmitter k sends its signal Sk to its intended 
receiver through beamforming vector f/^ G and that 

the Sk are statistically independent with zero mean and 
= 1 for all k G 1C. Under these conditions, the 
available baseband signal at receiver k before PS is ideally 
given by 

K 

Tk = hffcffcSfc -f ^ hgfjSj +nk, (1) 

desired signal ^ 

interference 

where hkj G denotes the channel vector between 

transmitter j and receiver k, and S C is the additive white 
Gaussian noise (AWGN) introduced by the receive antenna, 
which is assumed to have zero mean and variance a^. 

Each receiver splits its received signal to the information 
decoder and the energy harvester by means of a power splitter. 
Let pk (0 < pfe < 1) denote the power splitting (PS) ratio for 
receiver k, which means that a portion pk of the signal power 
is used for signal detection while the remaining portion 1 — pk 
is diverted to an energy harvester. Thus, the available signal 
for ID at receiver k can be expressed as 

4° = [ hffeffcSfc + -f nfc J Vk, (2) 

V j=t,jAk / 

where Vk is the additional AWGN circuit noise with zero 
mean and variance due to phase offset and non-linearities 
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during baseband conversion Q. We assume that each receiver 
employs single-user detection by treating the cochannel 
interference as noise. Then, the SINR at receiver k is given 
by 

T. Pkl^kk^kf 

k k - - > -- ~ T —■ 1-51 

Pk( E + ) +Wfc 

Besides, the total harvested energy that can be stored by 
receiver k is equal to 




ik (1 - Pk) 


K 




+ 


(4) 


\i=i / 

where G ( 0 , 1 ] denotes the energy conversion efficiency of 
the /cth EH unit, which means that a portion of the received 
RF signal is used for EH. 



Fig. 1. The X-user MISO interference channel. Each receiver .splits the 
received energy in two parts for ID and EH, respectively. 


B. Channel Error Model 

Because of many factors such as channel estimation 
errors, quantization errors, and feedback errors/delay, it is 
not possible in practice to obtain perfect CSl at both the 
transmitters and receivers. Let G denotes the 

estimated channel vector between transmitter j and receiver 
k. Then, the actual CSI can be expressed as 

hfcj = h.kj + Gfcj, J, k G K., (5) 

where e^j denotes the CSI error vector. To model the CSI 
errors, the well-known norm-bounded error (NBE) model 
is adopted, where we assume that the channel estimation 
error e^j is bounded in its Euclidean norm, that is 

\\ekj\\<Vkj, j,kGlC, (6) 

where rjkj is a known positive constant. Equivalently, h.kj 
belongs to the uncertainty set Ikkj defined as 

^kj = {h|h = hfej -f Bkj, llefcjll < 77fcj}. (7) 


The shape and size of model the kind of uncertainty in 
the estimated CSI, which is linked to the physical phenomenon 
producing the CSI errors. It should be emphasized that the 
actual errors {efcjj are assumed to be unknown while the 
corresponding upper bounds {rjkj} can be obtained using the 
preliminary knowledge of the type of imperfection and/or 
coarse knowledge of the channel type and its main charac¬ 
teristics 120). This model is particularly suitable for systems 


where CSI is corrupted by quantization |211. 


C. Optimization Problem 

We assume that the ID and EH units of each receiver are 
characterized by certain quality of service (QoS) and EH 
constraints. The QoS constraints require that the SINR of 
receiver k should be no smaller than a given positive target 
7 ^.. In the meantime, the EH constraints call for the harvested 
energy of receiver k to be no smaller than a positive threshold 
In this work, we focus on robust JBPS design under such 
constraints. Specifically, minimization of the total transmission 
power subject to the above two types of constraints, in the 
presence of imperfect CSI with NBE model, can be formulated 
as the following constrained optimization problem: 


min E llffcf 


Pk\ E 1+^2 




> Ik, 


( 8 ) 


s.t. 


K 


ik (1 - Pk) I^E | + crlj > ijjk, 
0 < Pfc < 1, llefcjjp < r)h, yj,kG 1C. 


Different from IIZ)> where it is assumed that the 

perfect knowledge of CSI is available, the goal of our work is 
to investigate the robust JBPS design, i.e. to guarantee that 
the SINR targets 7 ^ and EH requirement tpk are satisfied 
for all possible CSI errors. As compared to the non-robust 
design in 1T6),|T7|, the above robust JBPS design can provide 
guaranteed QoS and harvested energy for each receiver for 
all possible channel realizations in the bounded uncertainty 
regions Ikkj in (0. Solving the robust design problem ([^, 
however, is more challenging because there is an infinite 
number of constraints (due to the NBE model) and each of 
the SINR (EH) constraint is not convex. Both these properties 
make the problem (j^ very difficult to address. 

Similar to Lemma 3.1 & Lemma 3.2 in tn. it can be shown 
that the feasibility of ([^ is independent of the EH constraints 
and PS ratios. That is, problem (j^ is feasible as long as the 
following problem is feasible. 


find {ffe} 




> Ik, 


s.t. 


J = 1 1 J 7^ ^ 

llefcjf < pL-, yj,kG K. 


(9) 


This property provides an easier way to check the feasibility 
of (j^. In the rest of this paper, we assume (|^ is feasible. 


HI. Proposed Robust Design Based on Relaxation 


In this section, we propose two efficient approaches to 
solve the robust JBPS design problem ([^. In the first design 
approach, the celebrated SDR technique is applied to convert 
the semi-infinite constraints into linear matrix inequalities 
by means of the S-Procedure p2) , p^ . In general, the 
relaxed problem is not guaranteed to have a rank-one solution 
1161, p4) , p5| , we thus provide a good heuristic solution 
to recover a feasible solution to In the second design 
approach, we formulate problem ([^ into a SOCP problem 
based on SOCP relaxation. The resulting SOCP problem has 
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very low computational complexity and shows great potential 
for applications with large antenna arrays and large number 
of transmit-receive pairs. 

A. SDP with Rank Relaxation 


According to |26|, we introduce a new optimization 


K 


s-t- >al + ^ 


Ik 

K 




Pk 


E hf,F,h,, > 


ik{l-pk 


— a 


k ’ 


relaxed problem will be discussed in detail in Subsection III-B 


Pkj — 


max 


VeJ 




Qkj = mm 

Vef 




H 

kk 


-'fcfc 


'Ffc(h 


kk 


-kk 


^ Pkj+t^k + -^^ 


lefcfcf < vlk^k e /C, 


(h^- + e^)Fj (hj,^- + e^,j) < pkj. 


llefcjf < Pfcj, Vfc G lC,j ^ k. 

Similarly, the EH constraints in problem ( [T0| ) can be ex¬ 
pressed as 


(hfefc + efcfe)Ffc(hfcj, -f e;.^,) -f 


K 

j=td¥=k 


> 


i’k 

Cfc(i - Pk) 


2 


llsfefcf < ilk, yk G K, 


(hg 


-.H 


Fj (hfc,+e^ .) > qkj, 


\\ekj\r < vlp VA: G ICJ k. 


(16) 


variable F/j = ffcf^, V/c G /C and rewrite problem (|8|) as 
follows 

K 

min X) Tr(Ffc) 

{Ffc. Pk} k=l 


By applying the S-Procedure, the constraints in ( [T3| ) and 
111 can be reformulated to finite convex constraints, which 
are equivalent to (181 and shown at the top of the 

next page, where ak = and Afc^, Vj, fc G K, are slack 
variables. Similarly, we can recast and as and 


0 

Tnei 


with /3fc = 


-I ’ 

en, problem (|10|i can be expressed 


and slack variables ^j,k G 1C. 


( 10 ) 


{Fk 


s.t. 


min 

Qfci / 3fc, Afc, , Mfc, ,} • IPfe-i , QkA i^k 


as 

E Tr(Ffc) 

/c=l 


i=i 

0 < Pfe < 1, Ffc ^ 0, rank(Ffe) = 1, 

Ijefejf < T]lp Vj,fc G K.. 

This problem can be relaxed as a convex problem by dropping 
the non-convex rank-one constraint rank(Ffc) = 1, since the 
objective function and the constraints are linear in F^ and 
both convex function of pk- It is worth noting 
that the relaxation is not optimum and postprocessing of the 


([^ and 

ak > Y(3k > 1, inv^fc) + invp(/3fc) < 1, 
Pkj > 0, Qkj > 0, j ^ k, 

Fk h 0, Afcj > 0, Pkj > 0, Vj, keJC, 


(17) 


where constraints (18i-(21i appear at the top of the next 


However, this problem is still computationally intractable 
because it involves an infinite number of constraints. By 
applying the S-Procedure, the infinitely many constraints can 
be reformulated into finite convex constraints. 

We first observe that each term in the SINK and EH con¬ 
straints containing F^ in problem ( [T0| involves independent 
CSI errors. Hence, we introduce two auxiliary variables 


kGiC, j^k, (11) 

^kjFjhkj, keK, jy^k, (12) 


where pkj is the maximum (worst-case) cochannel interference 
power from transmitter j to receiver k and qkj denotes the 
minimum (worst-case) power available for EH from the 
transmitter j to receiver k. Then, with the help of these 
two variables, the SINK constraints in problem ( [T0| ) can be 
equivalently rewritten as 


(13) 


(14) 


page. In (17 1 , the set of constraints involving invp(-) must be 
satisfied with equality at optimality; otherwise the objective 
value can be further decreased by decreasing a^'s. The above 
problem is a convex SDP problem which can be solved by an 
off-the-shelf solver 113, ID. 

B. Proposed Rank-one Recovery Method 

The matrices Ffc obtained by solving the relaxed problem 
are not guaranteed to be of rank one in general, and 
hence, the solution provides a lower bound to the original 
problem (j^l- If F^; happens to be of rank one, then the 
principal eigenvector f^ of F/j, such that F^ = ffc ffc ^ will 
be the optimal solution to problem §□ Otherwise, one has 
to resort to other techniques to obtain a suboptimal rank-one 
solution from Ffc. In this work, inspired by we provide 
a simple heuristic, yet effective approach to overcome this 
difficulty when when higher-rank solutions are returned by 
solving problem (17 1 ^ 

Before we proceed to introduce the rank-one recovery 
method, we first calculate the worst-case channels for given 
beamforming vectors and PS ratios. Assuming that {f^} 
(the principal eigenvectors of {Ffc}) and {p^} have been 
determined in the previous subsection, then the worst-case 
CSI errors which minimize the SINK of user k, are the 
solutions to the following problems 

min |(hffc-feffc)f^| s.t. ||efcfcf < (22) 

{efcfc} 

max |(hg -f ekj)^*\ s.t. Ijefc^f < j ^ (23) 

i^kj} 


To solve the constrained optimization problem ( |22) l, we resort 
to the following inequality for any complex number x and y 


||a:| - IpII <\x + y\< ||a:| -f \y\\. 

According to ( [24| , we have the following inequality 

|hffcifel = |hffcffc+e«ffc|^ 

> llhffcffcl - leffcffcll > llhffcffcl -pfcfcllffclll, 


(24) 


(25) 


(15) 


f*|| = and fk is the largest eigenvalu e of Ffc. 


“^We can also use randomization techniques |26| to recover a rank one 
solution. 
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U 


ki^k: {Pkj}^kk: 


+ Afe/cl 


Ffch 


feilfc/c 


K 


^^kk^k ^h^^Fkhkk - E Pkj -(^l- ^lak - Xkkvlk 


7fc 


'^kj (Fj, Pkj , Xkj ) — 


A 


^^k{^kj \^Qkj\j^k^ ^kkT Pk^ 


— F + Afc,I 

H 
kj 

Ffc “ 1 “ pkk^ 


^1 ^ ^ J^'-kj 

~^kj^j Pkj ~ ^kj^j^kj — ^kjVk 


i=i j'v^fc 
—F ,h(i 

kj^j-i^j -■'‘j-ikj 

Ffchfefc 
K 


> 0 , j ^ fc 


h 0 (18) 


(19) 


^ick^k h^^Ffehfcfe + E Qkj ^^/l/c+O"^ PkkVkk 

j=i,J^k 


>- 0 


( 20 ) 


kj (Fj , Qkj } Pkj') 


A 


Fj “ 1 " Pkj^ Fjhfcj 

^kj Fj hg Fj hkj - Qkj - Pkj pIj 


h 0,j k 


( 21 ) 


where the second inequality follows from the boundedness of 
CSI errors and the Cauchy-Schwartz inequality. It is worth 
noting that we have used a mild assumption in (|25]) that 
Kkfkl > I which holds true in most cases since the 

CSI error is often much smaller than the channel coefficients 
and hkk cannot be orthogonal to ik- We observe that the 
first inequality in (p5]l holds with equality if and onl y if 


^kk^k — ^^kk^k’ ^ ^ 


L. Also the second inequality in (251 
holds with equality if and only if = /3f^, /3 G C. Thus 
the optimal solution e^.^, of ( 221 can be obtained as 

-Pkkh^k^k 




= conj 


I \^kk^k I 


(26) 


mKk^t 

^kk=m- (27) 

Proceeding in the same manner, the optimal solution of 


© can be obtained as 


a = 


Pkj 


H £•* 




= conj 


= /3f;. 




K^; 



(28) 


(29) 


*^kj — i-'^j ' 

Similarly, the worst-case CSI errors which minimize the 
EH of user k are the solution to the following problem 

min |(hf^.+eg)f*| s.t. ||efcjf < Vj, fc G/C. (30) 

X^kj f 


Let denote the optimal solution of (301 when j = k, 
which can be calculated following the same approach as 
used for (26 1 and (271. In the case j ^ k, we modify the 


where F* = f;(f;) 
associated with 


(31) 


cost function using the method of Lagrange multipliers 
which yields the following Lagrangian functiorj^ 

= {^kj + 6^)Fj(hj.j + ef.j) + Tkji^kj^kj ~ Vkj)^ 

j and Tkj is the Lagrange multiplier 
the bounded CSI error constraint from 
transmitter j to receiver k. Taking the gradient of C in ( [3T| ) 
with respect to conj(efc^), we can obtain 


-kj 


= -{F*+Tkjl)-^F*h,^ = 




(32) 

rkj + iitjir 

The Lagrange multiplier can be determined by solving 


^kj^kj 


Pkj = 0. 


(33) 


^In this case, |30| cannot be simply solved because the mild assumption 


does not hold in this case. 


Introducing gkj = becomes 


SkjSkj 


= Pk 


(34) 


Therefore, the lagrange multiplier Tkj can be given by 

= \/skjSkj/vlj - l|f;f ■ (35) 


We note that but in the case j 7 ^ k, and 

can not be simultaneously attained for the same channel 
realizations, which means that a CSI error vector minimizing 
both SINK and EH with given beamforming vectors and PS 
ratios does not exit in general. However, we employ both 
and ef,j to guarantee the robustness of the joint design. 

With the worst-case analysis described above, we recover 
the rank-one solution by scaling up the beamforming vector f^ 
by and then jointly optimize {</5fc} and receive PS ratios 
{pk} to satisfy both worst-case SINK and EH constraints and 
yet minimize the total transmission power. Specifically, we 
consider the following problem with given {f^} 


K 

min y]]<^fe||ffc 

{pfc, ‘Pfc} - 


/c=l 


S.t. 


Pk^k^kk 


K 


> Ik, 


E PkPjUkj + Pk<jl + ojI 

j-lj^k 


K 


£.k (1 - Pk) + ^kj > ^k, 

0 < Pfe < 1, tpfe > 0 , Vfc G 1C, 


(36a) 

(36b) 


(36c) 

(36d) 


where = ili^kj + denote the minimum 

interference power and = {|(hfj + eg)f;p} 

denote the maximum interference power. Let us define 
Xk = ^Ukk - Ylf=i,j^k‘P 3 ^kj - crl, Xk is implicitly 
supposed to be larger than 0 , and note that otherwise p 6 b[ ) 
will be infeasible. Thus, P 6 b[) can be rewritten as 


\\[2u}k,Xk - Pk]\\ <Xk + Pk- 


K 


(37) 


Similarly, by introducing yk = fj^kj + cr^, (36ci is 

i=i 

equivalent to 


\[2^/'^pk/^k,yk + Pk - l]|| <yk- Pk + 1- 


(38) 
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TABLE I 

Algorithm-1 : Proposed SDR-based robust design 


1. Solve problem (|l7| to obtain the optimal {F^} and {p^}- 

2. For each k, find the principal eigenvector {f^} of FJ. If 
rank(FJ) = 1, Wk G /C, then {f^} and {pj} are the 
optimum solution and exit the algorithm, otherwise go to Step 
3. 

3. Solve problem ||39[ with {f^} to obtain {(pk} {pk}- 

4. Return the beamtorming vectors {pfc}- 


Then problem ( [36) > can be reformulated as 


min 

{pk, pk} 


K 




S.t. I I Pfc]|| ^ Pkj 

\\[2y/llJk/^k,yk + Pk - i\\\ < Vk - Pk + 1, 

K ^ 

Xk = ^Ukk - E ^PjUkj - CTfc, 

3 = 1,j^k 
K 

yk= Y. VjUkj + <xli 0 < Pfe < 1, 

i=i 

ipk >0, Xk> 0, yk > 0, Vfc G K,. 


(39) 


The above optimization problem is a SOCP problem |30| 
because its objective function is linear and its constraints 
are linear or second-order cones. It can be efficiently solved 
by off-the-shelf algorithms. The proposed SDR-based robust 
algorithm for problem ^ is summarized in Table ^ 


C. SOCP Relaxation 


In Subsection A, we relaxed the robust JBPS problem as a 
SDP problem. It is well known that solving a SDP problem 
requires relatively high computational complexity. To obtain a 
low complexity solution, we propose to formulate the original 
problem as a SOCP problem based on proper relaxations, i.e. 
SOCP relaxation. Similar to ini’ problem (|^ can be relaxed 
as the following problem by replacing the EH constraints 
with the sum of the SINR and EH constraints 

K 


min 

E iiffcir 





s.t. 

(' + 4) 1 


K 




kj 


Pk ' 


(40) 


0 < p/c < 1, lle/cjjp < vlj, yj,ke K. 

According to (|2^, we have the following inequality 


|hg.f,| = |hgf,+egf,| 


< + \^kM < IlhfTjl +»7fejl|fflil. 3 ^ k. 


(41) 


Together with inequality ( |25| l, the SINR constraints in ( |40| ) 
can be reformulated as 


- PfefellffellP 

> E lihgf.l 

i=i 


■Pkjl 




(42) 


Similarly, the EH constraints in (40i can be expressed as 


1-f— 
7fc 


Ihffcffcl 


Pkk I ffc III — 


i’k 




Cfe(l-p/c) Pk 


(43) 


Eurthermore, by introducing auxiliary variables fikpi 
which satisfy 


hfe/cf/cl - Vkk llffcll > Pkk, 


|hfjf,l 


Vkj |l*-j 

(42 1 and (|43]l can be expressed as 


+ Vkj ||fj|| ^ Pkji j 7^ k, 


-Plk> 

Ik 


K 

j^lj^k 


2 I ^2 I 

+ <Xk + 




Pk 


1 + 


Ik 


Pkk ^ 


V’fc 


ikp-pk) Pk 


(44) 

(45) 


(46) 

(47) 


Let = Pk and h\ = 1 — Pk, it follows that a\ + b\ = 1. By 
further introducing > ^, and |h^fj| < ekj, 

we can write (|40|) as the following SOdP problem 


min t 

{ffc) ^k, ^k, (2k, dk, 0k j} , 

s.t. ||[ff,...,f^]|| <t, 

II [Jfc I Cfclll ^ 

Pfej||fj|| ^ + Pkj, |h^,jfj| < Bkj, j 7^ k, 

Pfefc||ffc|| < hffc ffc - I3k k, (48) 

ll[cfc, 4 ]ii < + ^Pkk,yi, 

\m'ii^ki^kf/\dk-bk]\\<dk + bk, 

II [2'ydUfc, Ck 0,k\ II ^ Ck 4“ U/c, 

ii[afc,(>fc]|| < 1, 

ak >0,bk> 0, Pkj > 0, Vj, k,elC 


where 


Jfe = [Pki, ■ ■ ■ fdk{k-i), /dk{k+i), • • ■, PkxY■ (49) 

The constraints ||[afc,5fc]|| < 1 must be satisfied with equality 
at optimality; otherwise the objective value can be further 
decreased by increasing a^'s. Note that we can restrict h^^ffc 
to be positive, which incurs no loss of optimality since we can 
always phase-rotate the vector f^ such that is positive 

real without affecting the cost function or the constraints. Since 
the solution to the SOCP relaxation problem ( |40| ) may not be 
a feasible solution to problem ([^, a robust solution recovery 
method must be employed to ensure the robustness, which will 
be introduced in the following subsection. 


D. Proposed Closed-form Robust Solution Recovery Method 

In Subsection B, we introduced a rank-one recovery method 
to recover a robust solution for problem ( [T7| by solving a 
SOCP problem ( |3^ . Due to the SOCP relaxation, the solution 
to problem ( |48| ) may not be robust to all possible channel 
realizations. In this subsection, we propose a closed-form 
recovery method in which {fj^} obtained by solving (481 are 
scaled up by a common factor y/p, and then jointly optimize 
(p and the receive PS ratios {pk} to satisfy both the worst-case 
SINR and EH constraints. Since we relax the EH constraints 
in (0, it is required that (p > 1 to ensure the robustness 
of the algorithm, which means that more power is needed 
to satisfy the worst-case EH constraints. It is worth noting 
that this closed-form method was first proposed in d) where 
perfect CSI is considered. 
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TABLE II 

Algorithm-2 : Proposed Robust Design with SOCP Relaxation 


1. Solve problem to obtain the optimal {f^} and {p^.}- 

2. Solve problem (32 with {f^} to obtain <p* and {pk}- 

3. Return the bearntorming vectors and {pk}- 


Similar to problem (36 1 , we consider the following problem 
with given {f^} 


min 

'p, {pk} 

s.t. - 

p 


K 

(fpkUkk 




> 7fc, 


(50) 


pkUkj-\-pk(7l-\-^l 

j = '^, 39 ^k 

Cfe(l - Pk) J2^Ukj + cr"^ > Pk, 

0 < Pfc < 1, > 1, \/k £ K., 

where the definit ions of Ukj and u^j have already been stated 

_ A ^ 

Introducing Xk = ^uuk - E '^kj and 

i=i 


in Section 


III-B 


K 

yfe = E ^kj, 

7 = 1 

rewritten as 


Vfc S JC, problem (50i can be equivalently 


mm 

V, {Pk} 

s.t. Pk > 
Pk> 




ipxk-crt 

4 


(51) 


ik{vyk+iyl) ’ 

0 < Pfc < 1, i-p > 1, yk £ 1C. 

Similar to problem (|5T| admits a closed-form solution, 
which is given by 


mm 

1 




s.t. ipk = 


1 , 


0<5fc(l)<l 

<p, gfe(l) > 1 orpfe(l) < 0 

P > Pk, V/c e 1C, 


(52) 


where gk{p) = 


Pk 


_ k I _ . 

(pXk-cr'i ik(<PVk+<x‘t,) 


and (fi is the largest root 


” -2- IS 

ip*Xk-a-l 


of the equation gk{p) = 1. Additionally, pk = 
the corresponding PS ratio with given p*, where p* is the 
optimal solution of ( [5^ . The proposed robust design with 
SOCP relaxation for problem ([^ is summarized in Table 
Remark: It is worth noting that the two algorithms intro¬ 
duced in Section III-B and III-D can both be employed to re¬ 
cover a feasible rank-one solution to the relaxed problem ( [T7| ) 
and ( [48| , respectively. The two algorithms exhibit a tradeoff 
between recovery accuracy and computational complexity. In 
Algorithm-1, we employ the rank-one recovery method from 
Section IIII-B while in Algorithm-2, the recovery method from 
is employed to further reduce the computational 


III-D 


Section 

complexity. It is also important to note that the two algorithms 
have very close performance in our simulations. 


IV. Proposed Iterative Robust Design 
Since the robust design approaches presented in Section 
l^rely on relaxation techniques, they may not guarantee the 
optimal solution, we propose an iterative algorithm based on 
CCCP ||3T|, @ in this section to improve the relaxation 
solutions and possibly obtain the optimal solutions. 


A. CCCP-based Iterative Robust Design 

The CCCP technique is widely adopted for solving non- 
convex problems | [3^ by transforming them into a sequence 
of convex programming problems. The basic idea behind the 
proposed CCCP-based algorithm is to iteratively approximate 
the original non-convex feasible set in ( [T0) i around the current 
solution by a convex subset and then solve the resulting convex 
approximation at each iteration p^ . 

While the conventional CCCP considers scalar functions, 
for the current application we need to extend it to the case 
of matrix functions. Our iterative method is motivated by the 
observation that the non-convex constraint rank(Ffc) = 1 is 
equivalent to the constraint 


ffeff = F 


k' 


(53) 


Hence, relaxing constraint (53 i to Ffe results in the 

following problem 


min P (54a) 

{Ufc, ffc , dfc, Pk, 2,kj , Pkj}, {Pkj , Qkj} j,/,k , t 

S.t. ||[ff,...,f^]|| < f, (54b) 

([^, ([2^ and ([2^, (54c) 

ttfe > 1, > 1, invp(afe) -f invp(/3fc) < 1, (54d) 

Pkj > 0, Qk] > 0, jV k, (54e) 

ffcf/f ^ Ffc, (54f) 

Ffc ^ 0, Xkj > 0, Pkj > 0, Vj, k £ 1C. (54g) 


Lemma 1: Problem ( [T0| ) is equivalent to problem (54i. 
Proof: As we can see, problem ( [54| is a relaxed version 
“ :,^of 

= Ff, we can assert that problem (lOi is 
equivalent to problem (54 1 . We prove f^f^^ = F^ by showing 
that rank(F^) = 1, which is quite obvious since fjJfE 


of problem (101, thus if the optimal solution {f^ 


satisfies 




rank one matrix. This completes the proof. 

As we can see, only constraint ( |54f| ) in problem ( [54| i is 
non-convex. Thus, with the help of the following inequality 


ffcff - (f^(f, - + (ffc - f^)f^^ - 

= (ffc-f^)(ffc-f^)^^0, 

we can transform ( |54f| ) to 

^f^(ffc-f^)^ + (ffc-f^)f^^ + f^ 


ffcff 




h Ffc, 


(55) 


(56) 


where denotes the current feasible point in the ith iteration. 

With ( [56] ), problem ( |54) i can be reformulated as the follow¬ 
ing convex optimization problem in the ith iteration of the 
proposed CCCP-based algorithm 


min P{^k) 

{Ffc. f;,. ^k^ 0k j ^kj^ P'kj}t {Pkj^ Qkj}j^k^ ^ 


s.t. II [f' 


T 


> fA:]ll ^ 


(^,(1^,(1^ and (1^, 


invp(afc) -I- invp(^fc) < 1, 

Pkj > 0, Qkj > 0, jV k, 

f^(ffc - fl)^ + (ffc - + f^f^" h Ffc, 

Ffc ^ 0, Xkj > 0, Pkj > 0, Vj, k£lC, 


(57a) 

(57b) 

(57c) 

(57d) 

(57e) 

(57f) 

(57g) 

(57h) 


































TABLE III 

Algorithm-3 : Proposed CCCP-based Iterative Robust Design 


Define the tolerance of accuracy (5 and the maximum iteration 
number Initialize the algorithm with a feasible point 

{f^, pfc}- Set the iteration number i = 0. 

Repeat: 

- Solve problem i |57| with the current feasible point 

- Assign the solution to {f^"^ , p/c} update the iteration 
number : i = i 1. 

Until: the objective function converges, i.e. — P*| <5 

or the maximum number of iterations is reached, i.e. i > 


where P{fl) = 

Lemma 2: Suppose problem ( [57| is feasible, then strong 
duality holds true for problem (|57|l and its dual problem. 

The proof is relegated to Appendix 

To summarize, we can see that the feasible set of problem 
([54|i is a subset of the original set defined in problem ([T7|. 


Then, if the initial point {f°} is feasible for (54i, all the 


feasible points obtained by iteratively solving problem 
( [54| i always belong to the true feasible set of ( [T7] i. The 
proposed CCCP-based iterative robust design for problem ([^ 

Regarding its convergence. 


is summarized in Table III 


yPfcl(hffc+effc)ffcl 


VPl\ +<yl ]+uil 


> 


Pfcl(hfcfc+efcfc)ffc I 


Pl\ +(tI ]+uil 


> 1 . 


Let 


K 


= 6(1 - Pi) 






-'kj 

expressed as 


where ejf, is calculated by (|32|| and (|35|l. Then (p can be 


P = 


li’h 


- ut 


(60) 


f;.(i-p^) 

We remark that Algorithm-3 consists of a two-stage algo¬ 
rithm for solving problem (|^. In the first stage, the initializa¬ 
tion method is applied to find a feasible solution of problem 
(H). If the initialization method fails to find a feasible solution, 
Algorithm-3 declares failure and stops. In the second stage, we 


iteratively solve problem (57i. 


we 


have the following Lemma. 

Lemma 3: Algorithm-3 produces a non-increasing sequence 
of objective values. Moreover, every limit point {f^} of the 
iterates generated by Algorithm 3 is a KKT point of problem 

Proof : Please see Appendix 

B. The Proposed Initialization Method 

The proposed Algorithm-3 presented in the previous sub¬ 
section requires an initial feasible point of problem ® Eg. 
If problem (57 1 is initialized with an infeasible point, then the 
CCCP may fail at the first iteration. 

The proposed initialization method is based on Algorithm-2. 
As has been stated in Section |III-C| the solution obtained by 
solving (48 1 may not be robust to all channel realizations due 


V. Complexity Analysis 

In this section, we compare the relative computational 
complexities of the proposed robust design algorithms. As will 
be seen from our analysis and the simulation results in Section 
VI the proposed robust algorithms exhibit a tradeoff between 
computational efficiency and robust performance. For ease 
of comparison, we assume that all transmitters are equipped 
with the same number of antennas, i.e., Nk = N, Vk G 1C. 
Moreover, we apply the basic elements of complexity analysis 
as used in EH- 

1) Algorithm-1: Consider problem ( [l7] i, which involves 2K^ 
linear matrix inequality (LMI) constraints of size N -\- 1 
and K LMI constraints of size Here, the number of 
decision variables n is on the order of 0{KN'^ G-AK"^). Thus, 
the complexity of a generic interior-point method (IPM) for 


solving problem (17 1 is on the order of the quantity shown 


to SOCP relaxation. In order to make the initialization method 
as simple as possible, we propose a new method instead of 
the robust solution recovery method in Section |III-D| We 
first observe that the solution obtained by solving^T^ can 
provide guaranteed SINR levels if ( [48] l is feasible. Thus, if 
we scale up the beamforming vector {f^} by a sufficiently 
large common factor then will be a feasible 

point of problem (|^, since 


on the first row of Table IV Regarding the complexity of 
the rank-one recovery method ( (39] l, we note that it involves 
2K variables and 2K second-order cone (SOC) constraints 
of dimension 3. Hence, it follows that the complexity of the 
SOCP problem is 0{2KV^ [2K3^ + AK^]). We note that 
the complexity of the worst-case channel vector calculation 
is dominated by singular value decomposition and matrix 
inversion operations, and it is considered to be negligible 
compared to solving problem |T7|. 

2) Algorithm-2: Problem (48 1 involves 2K^ -\- AK -\- 1 
SOC constraints, including 1 SOCs of dimension KN -\- 1, 
K SOCs of dimension K + 2, SOCs of dimension 
N -\- 1, AK SOCs of dimension 3, and — K SOCs of 
dimension 2. The number of variables n is on the order of 
0{KN -\- 3K A- 2K^). Thus, the complexity of a generic 
IPM for solving problem ( [48] ) is the quantity shown on the 
second row of Table IV Here, the complexity of the closed- 


form robust solution recovery method is dominated by the 
worst-case channel vector calculation and solving a quadratic 
equation, which is negligible. 


(58) 


(59) 


3) Algorithm-3: Problem (57 1 involves KN"^ -\-KN -GAK"^ -\- 
2K -\- 1 variables, 1 SOC constraints of size KN -\- 1, 2K^ 
LMI constraints of size W -f 1, and 2K LMI constraints of 
size N. Hence, the complexity of the CCCP-based iterative 
robust design is as shown on the third row of Table 

We can check the asymptotic complexity of the proposed 
algorithms when N and K are large, i.e. let V = AT — oo. 

^Here, we ignore the constraints of lower sizes, since they will not affect 
the order of the whole problem. 
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TABLE IV 

Complexity analysis of the robust designs 


Robust design 

Complexity Order (suppressing the ln(l/£:) ) 

Algorithm-1 

0{ny/2K^(N -1- 1) + KN [2K'^{N + 1)'’ -b KN^ + 2nK^{N -b 1)-^ -b nKN^ -b n^]) 

■yO{2KVtK[2K3'^+4:K'^]), n = 0{KN'^ + AK^) 

Algorithm-2 

0{ny/AK'^ +%K + 2 [{KN + 1)^ -b K{K + 2)^ + K'^{N + l)^ + AK2A + (K^ - K)2'^ + n^]) 

n = 0{KN + 3K A-2K'^) 

Algorithm-3 

The complexity of Algorithm-2 + 0{nN'^^^■^y2K‘^{N + 1) + 2KN + 2 [2A'^(A^ + 1)^ + 2KN^ 

-\-2nK^{N + 1)2 + nKN^ + (KN + 1)^ + n = 0(KN^ + KN + 4K^ + 2K-\- 1), is the iteration number 


One can verify that the complexities of the proposed algo- 
are on the orders of 66A^^ and 


IV 


rithms in Table 

respectively. As seen, the lowest complexity 
is achieved by Algorithm-2, followed by Algorithm-1 and 
Algorithm-3. 


VI. Simulation Results 

In this section, we evaluate the performance of the proposed 
robust JBPS algorithms numerically. We assume there are K = 
3 transmit-receive pairs and all transmitters are equipped with 
Nk = N, fc S {1,2,3} antennas unless otherwise specified. 
We assume that each transmit-receive pair has the same set 
of parameters, i.e., 7 ^ = 'j, tpk = Cfc = 

= 09 ^ and rjkj = rj, Vj, k G K. for simplicity. Moreover, the 
nominal channel vectors {hfej} are randomly generated from 
independent and identical Rayleigh fading distribution with 
average power 1. We set ^ = 1, = —30 dBm, = —20 

dBm, (5 = 10“^ and = 20 in all our simulations. All 

the modelling and solution of the algorithms are performed 
using CVX p5| on a desktop Intel (13-2100) CPU running at 
3.1GHz and 4GB RAM. 

1) Feasibility rate: We first present the feasibility rates of 
the three robust JBPS design algorithms. In the simulation, a 
robust design algorithm is considered infeasible for a channel 
realization if CVX reports an infeasible status or 5;^ < 0 in 
the robust solution recovery method. The feasibility of the 
non-robust design are tested with 100 channel errors 
satisfying the NBE model for each channel realization. Fig. 
I^and Fig. [^present the simulation results obtained over 1000 
channel realizations. One can observe from this figure that the 
three algorithms exhibit similar (almost identical) feasibility 
rate compared to the boundj^The non-robust method fails to 
satisfy both the SINR and FH constraints almost all the time 
under NBF model. 

2) Transmission power: We illustrate the performance of 
the three robust designs in terms of average transmission 
power over 1000 problem instances. Fig.|^shows performance 
comparison among the three robust designs, the transmission 
power being averaged over problem instances where the robust 
designs are all feasible. It is observed that, as a price paid 
for guaranteed worst-case performance, the robust designs 
require higher average transmission power than the non- 
robust design. However, Algorithm-1/2/3 show near-optimal 
performance compared to the bound. The performance of 
Algorithm-3 is slightly better than Algorithm-1 due to the fact 
that higher-rank solutions may be returned by solving problem 

^The feasibility rate of the bound is equivalent to the feasibility rate of 
solving problem mi without rank recovery. 



Fig. 2. Feasibility rate (%) versus various 77 . = 4, 7 = 10 dB, '0 = 5 

dBm. 



Fig. 3. Feasibility rate (%) versus various N = 4, 77 = 0.1, 0 = 5 dBm. 



Fig. 4. Transmission power versus SINR tai'get 7. A/" = 4, 77 = 0.1, 0 = 5 
dBm. 
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Fig. 5. Average worst user SINR versus SINR target ■y. N = 4, rj = 0.1, 
'ip = b dBm. 



Fig. 7. Average worst user harvested energy versus EH target N = 4, 
T] = 0.1, 7 = 10 dB. 


0- In the low SINR region, almost all the solutions returned 
by Algorithm-1 are rank-one (optimal), thus the performance 
of Algorithm-3 is slightly inferior to that of Algorithm-1 since 
Algorithm-3 does not converge to the optimal solution in 
j^max iterations. 

The average worst user SINR performance of the proposed 
robust designs is illustrated in Fig. for both the robust and 
non-robust designs under the same simulation parameters as 
in Fig. 1^ Clearly, the average achieved worst user SINR of 
the robust designs are all above the SINR target while the 
non-robust design fails to satisfy the SINR constraints. 

Fig. § shows performance comparison among the robust 
designs for various EH constraints. One can see that the 
robust designs require higher average transmission power than 
the non-robust design. The best performance is achieved by 
Algorithm-3, followed by Algorithm-1 and Algorithm-2. The 
average worst user harvested power of the proposed robust 
designs is illustrated in Fig. where the simulation parameters 
are the same as in Fig. Clearly, the average achieved worst 
user harvested power of the robust designs are all above the 
EH target while the non-robust design fails to satisfy the EH 
constraints. 



Fig. 6 . Transmission power versus EH target Tp. N = 4, r] = 0.1, 7 = 10 
dB. 

3 ) Execution time: We then compare the performance of the 
robust and non-robust designs in terms of average execution 
time over 20 channel realizations. Eig. demonstrates the 
execution time (on a logarithm scale) versus the number of 


transmit-receive pairs with fixed number of transmit antennas 
N = 18. It is observed that the time consumed by all four 
algorithms increases with K. However Algorithm-2 consumes 
much less time than Algorithm-1 and Algorithm-3, which 
means that Algorithm-2 shows great potential for applications 
with large antenna arrays and large number of transmit-receive 
pairs. Algorithm-3 requires the most time as a price for better 
performance when higher-rank solutions are returned. 



Fig. 8 . Comparison of execution time versus K between the robust and 
non-robust designs with fixed N = 18. 'ip = b dBm, rj = 0.1, 7 = 10 dB. 

VII. Conclusion 

In this paper, we considered the robust JBPS design problem 
for multiuser MISO interference channel under a NBE model 
for the CSI. Three different robust design approaches were 
proposed to handle the highly non-convex JBPS problem 
with different performance and complexity. In the hrst design 
approach, we proposed to relax the original problem as a SDP 
problem based on SDR, which provides a lower bound for the 
robust JBPS problem if a rank-one solution is returned. A rank- 
one recovery method was provided to obtain a feasible rank- 
one solution if a high-rank solution is returned. In the second 
design approach, we formulated the robust JBPS problem as 
a SOCP problem based on SOCP relaxation and the Cauchy- 
Schwarz inequality. As compared to the SDR-based algorithm, 
the SOCP method has lower computational complexity, while 
achieving a performance very close to the performance bound. 
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We also provided a closed-form robust feasible solution re¬ 
covery method. In the third design approach, a CCCP-based 
iterative algorithm was presented to achieve near-optimal 
performance when a higher-rank solution is returned by the 
SDR-based algorithm. It was proved that any limit point of 
the iterative algorithm is a KKT solution to the robust JBPS 
problem. The simulation results showed that the proposed 
robust transceiver designs have near-optimal performance in 
the presence of imperfect CSI. 

Appendix A 
The Proof of Lemma[3 

We introduce the following variables to dualize the cor¬ 
responding constraints in problem ( |57] l (the notation A : B 
denotes the constraint B and its dual variable A) 


Ufc = 


Xfc = 


UP 


H 


Ufe Ufc 


Uk 

Xfe 

Xk 


xp 


H 


: (181, Ykj = 
: (20), Ykj = 


■ ^kj ^ O 7 ^kj ■ l^kj ^ O 7 
Ckj ■ Pkj ^ O 7 dfi^j . Qjij ^0, j ^ 


^ kj V ]^j 

Ykj Ykj 
Ykj Vki 


: (19), 

: ( 21 ), 


e : (576), ak : Ofc > Pk ■ /3k > 1, hk : (|57e|), 

Zk^ho, Ak-. (ra. 



akk = Ukillk - Tr(Ufc) 

_ -p f Tfc 

“ ai||hfefeh« II lai’ 


'ykVkk iTrlTur 1 

2aiArfc||hfcfch» ||1 '' 

IkvikNk _ Ikvlk 


ai||hfcfch(;(;,|| 2aiA''fe||hfefeh(;(,|| 2ai ||hfefeh((j, | 


> 0 , 


A 


akj = Vkjvlj - Tr(Vfej) 

— mini — , 

^ no ■ 


^kj 


^kj 


> 


aallhfcjh" I 


a2 ’ 2a2Affe||hfc>" 
■nljNk _ nl- 


a2||hfcjh((j. II 2o2Arfc||hfc3h(; 


}Tr(I^J 
> 0 . 


A 

^kj — ^k '^kj — 


> 


7fc 


ai hfcfch" 


ai||hfcfchp^|| 

7fc 


= 0 , 


A 


^kj — ykj Xk 
1 


ai(if-l)||h,,h"| 


> 


a3(Ff-tl)||hfcfch" II 

^ > 0 , 


ai(if-l)||hfc,h"|| ai(K-ri)||hfc,h« 


~ A 2 a ^ 

Oik — — X}^-^ , 


(69) 


and Zfc can be expressed as (70 1 . which implies that {U/j}, 
{Vfcj}, {Xfe} and {Y/jj} define a strictly feasible point of 
the dual problem and the resulting dual problem is a bounded 
problem. Thus, the dual problem is always strictly feasible. 
Together with the fact that problem ( |57) i is feasible, we can 
see that Slater’s condition always holds for the dual problem 
( |62| ). Hence, problem ( |57j l can attain its minimum and strong 
duality holds p2). This completes the proof. 


Appendix B 

The proof of Lemma[3] 

In this appendix, we extend the convergence proof of 
to the case of matrix function. We denote problem ( [57| ) by 
VXk). Define R{ik,Fk) = ffcff - and R(ffc,f^,Ffc) = 
tlifk - (iP + (ffc - - Ffc. It follows that 

R(ffc, ffc, Ffc) = R(ffe, Ffe). In the following, we complete the 
proof through three steps. 

In the first step, we show that each {ffc},i = 1,2,... 


is feasible to problem (54i. To this end, it suffices to show 
that is a feasible solution of problem (54l, provided 


(61) 

Then the dual problem of problem (57 1 can be expressed as 

max inf C (62) 

r z.to,u.to.v.,to,x.bo. 

< 'YkjhO,akj>0,bkj>0,Ckj>0,dkj>0 
[ Zfc^0,e>0.3fc>0./3fe>0,/ifc>0 

where C can be expressed as ( |6^, shown at the top of the next 
page. A^ume ||hfefch^j,|| =Jikj\\hkj'^\\, hkj > 0, then if we 
choose Ufc, Vfcj, Xfc and Ykj as in (64i, V/c G JC, j k, with 
02 > Oi^, 03 > and 1 > ^ + 01 , 02,03 e 

K+, and choose Afc = I, e and hk as any positive real number, 
we have 


that {f^} is feasible. Assume {f^} is feasible to problem (54i, 
thus we have R(f*., f*, Ffc) = R(f*, Ffc) = f^f,^^-Ffc A 0. It 
follows that there must exist {ffc^^} that is feasible to problem 


7^(f^), that is, {f, 
follows that 


'i+i 


} is such that R(fj 


i+l 


f^Ffc) ^ 0. It 


R(f*+\Ffc) = - Ffc h R(f^+\f^,Ffc) 


kJ — 


(fk-m, 


iH 

k 


Pkfp -FkhO, 


(65) 


( 66 ) 


(71) 

where the first matrix inequality comes from ( [55] ). This implies 
that {ffc^^} is feasible to problem (54i. This completes the first 
step. 

In the second step, we show that the objective value 
sequence {F(f^.)} monotonically decreases as the iteration 
index i increases. We denote the optimal solution to 7^(ffc) 
in the ith iteration by {f^,F^,p®}- According to the hrst 
step, {f^.,F^,p®} is a feasible solution to V{fl). Moreover, in 
the {i -f l)th iteration, the solution {f^+^, F^+^, is the 

optimal solution to 7^(ffc). Thus, we have < P(ffc), 

implying the monotonic convergence of {F(f^)} since P(ffc) 
is bounded below and P(ffc) is convergent. This completes the 
second step. 

In the third step, we prove that any limit point {f^, Ft, pti 


of the iterates {ft,Ft,pt} ^ KKT point of problem (lOl. 

Let S(ft) and C(ft) denote the solution set and constraint 
set of problem vlfl). We first prove e S(ffc)- Since 


Following the same derivation, we can obtain bkk > 0 and 
bkj > 0. Furthermore, we have 


(67) 


( 68 ) 


{ft! is a limit point of {f^}, there must exist a conver¬ 
gent subsequence {ik} such that limj_>ooffc^ = ffc- Since 
the objective function P(ffc) of V{fl) is strictly convex in 
{ffc S the point {ffc^^} is unique Hence, the 

entries of the two sequences, {F’(f^)} and (f^, have a one- 
to-one correspondence. By restricting to a subsequence, we 
can assume that converges to a limit point f^*. 

Dehne the constraint set C^(ffc) = {ffc|R(ffc, ffc, Ffc) 0} 
andjC^(ffc) = {ffc|R(ffc,ffc,Ffc) ^ 0}. It follows that 
C^(ffc) C C^(ffc), Vffc. Let us consider the set (ffc*)- Since 
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C = t^+f: Tr{[ Afc - iUfc + E - Xfc - £ Y,fc - ]F4 + £ 

A;=l j^k j^k k—1 

K _ K K _ K K 

+ J2 (-Ofcfe - Tr (Ufc) + Ukvlk)^kk + J2 J2 {-o-kj - Tr {Ykj) + VkjVkj)^kj + J2 J2 (-^kj + Uk - Vkj)Pk] 

k—1 k—lj^k k—lj^k 

K _ K K _ K K 

H“ ( ^kk Tr -\- Xk'i^i^f^)l-i-kk “1“ ( ^kj Tr j -h ykjTj^AiJykj ( dkj Xk Ukj^Qkj^ 

k—1 k—lj^k k—lj^k 

^ 9 ~ t ^ ~\ ~ ~ 

+ x; (ufcWfc - ak)ak + [Xk^- /3fcj Pk + ak + Pk + hk{mvp{ak) + invp(/3fe) - 1) 


K 


+ ^Tr{Ak{il{ik - flf + (ffe - - e Tr(ffeff) - tj , 

Ufc = Ufe + Ufchf^ + hfc^uf + Vjfe = Vjfc + + h^k^fk + Vjkh^k^fk^ 

Xfe = Xfc + xfchffc + hj^fcxf + Yjfc = Yjk + Yjkhfk + ^jkVfk + yifch^fchj^- 


K 


Ufe = 
Vfc,= 

X. ^ 

Y. ,= 


Ufc Ufc 


U 


H 


_fc '^k 
^ kj V kj 
_^kj vkj 
Xfe Xfc 
X_f Xfc 
Yfcj Yfej 

L ykj ykj j 


IkTlkk 






=‘|“n"i"iA-2„».i!h„K«i 


^ a 


7fc 


= diag < min 


^kk 


I |^03(i<'+l)’ 2a3Affc(-^^+l)||hfefch" II j ’ 03 (if+ 1 ) || 


= diag 


iag|min|^ 


^kj 


i||L 

t.hffc j 


1 

02 ! 

II 

1 


.11 J 

a,. 


1 1 


>0, 

>- 0 , 


ai(K-l)\\h.kjh.l 


^ 0 , 
^ 0 , 


(64) 


=I-^Ufe 

Yt. ^ 


Ik 

1 


I - 4-Ufc - Xfe - V Y,fe = I - 4- mini3^,-4^2*^in 

7fc ■5'' 7k l-ai ’ 2aiArfchfcfch« J ^'''= 


£ V,fc - Xfc - £ Y,fc h 

j^k j^k 

a3(-fS'+l) ’ 2a3Nk(K+l)hkkh^k ’ 2ai(if-l)Affchjfch» 

- 1 ^ - o/, h,, vi^. . ^ 1 _ 


K 


hfcfchf^. _ ^ 

, = r1^fc ai(if-l)h,fchffc 


_ 

Cff J *-K ^ 

j=l:3^k , , . 

V J_LfTfeJi Tfchfcfchl^'i _ / I t i_ 

7;=''ai aihfcfchfE ''“3(-fS'+l) a3(if+l)hfcfeh» 

4 n _ 


iitk 

'^' ik _ hfcfchU 

NkhjkhfJ a3(if+l)hfcfch« 


. ^ 0*1 a3(K+l))I^0 

3 = 1,O^k Ik 


(70) 


R(ffc,ffc,Ffc) is continuous in and lim^-^ooffc^ = then 
there must exist, for any fixed G C^(fjJ), an integer /f^ 
such that 

R(ffc,f;^Ffc)^0,VJ>n. (72) 

This implies that there must exist a sufficiently large T such 
that 

(C) C (f^o C (f^o , Vj > T. (73) 

Since G S(ffcO’ we can see that 

P{fk) > P{fp^"), Vffc G C^(ffc*) C C^(f;7. (74) 

Moreover, since P( ) is a continuous function, we have by 
letting j —^ oo in that 

P(ffc)>P(fr)i VffcGn(ffc*). (75) 

It follows from the continuity of R(ffc, f/^, F^.) that 

7^(ffc)>P(fr), VffcGn(ffc*). (76) 

One the other hand, we can infer from the second step that 

p(4*) = p(fr). (77) 


problem (54 1 and 
C^(47- h follows 
we 


Furthermore, since 4^ is feasible to 
R(4^4^Ffc) = R{fp), we have fp G 
that 4 G C^(4). Combining this with (|76| and (f77]l, 
obtain 4 G S"(4). 

According to Lemma we have 4^*^ = Then we 
argue that (4, F^, p^} satisfy the KKT conditions of problem 
(54i. With Lemma 1^ and 4 S §(4)’ there must exist optimal 
Lagrange multipliers as follows (fix ak,Pk) 


(78) 


A 

[ u; 

K 

Vf ■ = 

> kj 

[v- 

ylf 

* 

^kj 

'^kj - 

A 

< 

A ■ 

x*k . 

■V"* — 

’ ^ kj — 

y*k3 

y*k3 J 


■^/c’ 


a 


kji ^kji ^kj") ^kj‘ 


Together with {F*, f*, a*,/3*, A^-,g^-}, we can si 
that the following KKT conditions of problem 7^(4) hold 


Ufe(F*,{p4}^.^^,A;fc,a*)^0, 

Vfe,(F*,p4,AL) PO, k, 

'^kjXjjQkj’ kkj) 

F^ ^ 0, A4 > 0, > 0, VL G /C, 


( 79 ) 
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I - +<h,,hf,)+ 

E (V*fe + ^jk^fk + hjfeV^f + ?;*fchjfchffe)- 

i/fc ^ ^ ^ ^ 

(Xj, + + a;feh^^,h^) — 

E (Y*fc + y*fchffc + ^jkYjk + ylfehjfchfJ ^ 0, 

i/fc _ _ _ 

U* h 0, V* h 0, X* ^ 0, Y* ^ 0, 

<k = <4k-Mvl)>0, 

a*kj = vljvlj - Ti-(Vfc^) > 0, j ^ k, 

b*kk=xlvlk-^<X)>0. 

Kj = vljVlj - Tr(Y;^.) > 0, k, 

4, = > 0, 4, = - ^*k > 0, 

Tr(U*U,(F*,K^.}^.^^, A*,„a*)) = 0, 

Tr(V^^.Vfe,(F*,p*fc^.,A^,)) = 0, j ^ k, 

Tr{XlXk{Fl,{q*k,}^^,,^^lk,P*k))=0, 

Tr(Y*^.Y,,(F*,q*^,A.*^.)) = 0, J ^ k, 

Tr(Z*F*) = 0, 

aLK, = 0 > b*k^h = 0 . 

4jPlj=0> dl^qt^=0, j^k, Vj,fce/C, 


(80) 


( 81 ) 


where we have used the fact = F^. Note that (791 

denotes the primal feasible conditions, (80 1 denotes the first- 
order necessary optimality conditions and dual feasibility 
conditions, ( (8T] l denotes the complementarity conditions. Eqs. 
(79i-(81 1 imply that {f^,F^,p^} is a KKT point of problem 
(54 1 . From the above analysis, we conclude that {f^,F^,p^} 
is a KKT point of problem ([T0|). This completes the proof. 
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